EMPLOYING COMBINATION PROCEDURES TO 
SHORT-TIME EOP PREDICTION 
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ABSTRACT. A well known problem with Earth Orientation Parameters (EOP) prediction 
is that a prediction strategy proved to be the best for some testing time span and prediction 
length may not remain the same for other time intervals. In this paper, we consider possible 
strategies to combine EOP predictions computed using different analysis techniques to obtain 
a final prediction with the best accuracy corresponding to the smallest prediction error of 
input predictions. It was found that this approach is most efficient for ultra-short-term EOP 
forecast. 
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■ 1. INTRODUCTION 
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Prediction of the Earth Orientation Parameters (EOP) is a practically very important and 
theoretically very interesting task, one of the main fields of activity of operational EOP 
services. Various methods are developed to compute a highly accurate EOP forecast. How- 
ever, a usual and well known problem is that different methods show different accuracy at 
different time intervals and prediction lengths. A method which is the best for short-time 
prediction may not be such for long-time prediction, and vice versa. On the other hand, a 
method which was proven to be the best for a testing period of time may not remain the 
same for the currently considered period. As an example, results of actual Polar Motion 
(PM) and Universal Time (UT1) predictions are shown in Figure [TJ In these plots, three 
predictions compared are made at U.S. Naval Observatory (BA — IERS Bulletin A), at 
Pulkovo Observatory (AM — prediction of the IERS Bulletin A data), and SI — prediction 
of SLR series computed an the Institute of Applied Astronomy^. 

From this example, one can see that different predictions show the best results for differ- 
ent years and lengths of prediction. This is true not only for different prediction methods, 
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1 ftp://quasar. ipa.nw.ru/pub/EOS/IAA/erp_rs.dat 
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Figure 1: Results of actual UT1 predictions for three years 2003, 2004 and 2005. BA - 
Bulletin A prediction, AM — prediction of the same NEOS series made by the author, SI - 
author's prediction of the SLR series computed at the Institute of Applied Astronomy. 

but also for variants of the same method. All the methods practically used for EOP predic- 
tions use various configuration parameters, such as the length of reference interval, the order 
of polynomial and number of harmonics for least squares adjustment, or the orders of au- 
toregression and moving average for ARIMA. A prediction procedure adjusted for some time 
period for computation of the best EOP forecast may not remain the best for the following 
time span. Evidently, the main reason of this is that the Earth rotation is more complicated 
process than we are able to describe by our forecast models. 

So, the question is, whether we can develop a prediction strategy robust to irregular 
behavior of the Earth? One of the possible ways to achieve this goal is making use of 
combination procedures applied to a set of predictions made with different methods. Such 
an approach was tested during processing of the results of the EOP prediction comparison 
campaign (Luzum et. al, 2007; Schuh et al., 2008, Kalarus et al., 2010), and showed a 
promising result. Our purpose was to perform further investigation of possibility of improving 
the accuracy of EOP prediction using various combination procedures. In this paper, we 
have performed several tests of computation of combined predictions of PM and UT1 based 
on a posteriori and actual predictions. 

2. TEST DESCRIPTION AND RESULTS 

Three combination strategies proposed in (Malkin, 2009) were examined in this study using 
different input data: 

CI — the average of input predictions as proposed and tested by Luzum et al. (2007), and 
Schuh et al. (2008). 

C2 — the best previous prediction for given length. Normally, we compute EOP prediction 
on the day of the last observed epoch ( "today" ) . To compute n-day ahead prediction 
we examine the set of predictions made n days ago and select one that predicts today 



EOP most accurately. Then we compute our today n-day ahead prediction making 
use of the method used for computation of the best prediction made n days ago. 

C3 — the best yesterday prediction. To compute this prediction we examine the set of 
predictions made yesterday and use for today prediction the method which corresponds 
to the prediction made yesterday which predicts today EOP most accurately. 

Note that C2 and C3 predictions also use CI results along with the input data. Strictly 
speaking, C2 and C3 are not really combinations, just selections of the method that assures 
the best previous prediction is used to compute the current one. These three combination 
strategies were tested using two kinds of input data. 

The first set of input predictions consists of nine daily prediction series computed from 
the IERS C04 final EOP series starting from shifted epochs from 1 October 2006 through 
31 December 2007. Predictions were computed making use of prediction strategy developed 
at the Institute of Applied Astronomy (Malkin & Skurikhina, 1996). The method is based 
on use of three prediction techniques: 

1. Least squares fitting for trend and several harmonics (LS); 

2. Autoregression (AR); 

3. Autoregressive Integrated Moving Average (ARIMA). 

Final prediction is composed of several segments for different prediction length intervals 
computed using different methods and/or model parameters and merged in one continuous 
series. LS+AR combination was found to be the most effective for UT1 prediction, and 
LS+ARIMA combination proved to be the best for PM prediction. Later, such an approach 
was developed by Niedzielski & Kosek (2008) The main model parameters that can be varied 
to adjust the prediction procedures are the following: 

1. base interval used for LS fitting; 

2. number and periods of harmonics; 

3. trend order; 

4. AR order; 

5. IMA (Integrated Moving Average) order. 

Thus, nine prediction series were computed using various sets of listed parameters aside 
from item 2 which remains the same for all the predictions. Then three combined predictions 
CI, C2, and C3 were computed for the period from 1 January 2007 through 31 December 
2007. Figure [2] shows the result of comparison of input and combined predictions with the 
IERS C04 EOP series. Both the RMS and maximum differences between them are depicted 
in the plots. 
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Figure 2: The RMS and maximum errors of the input and combined predictions: a posteriori 
input predictions. 
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Hereafter we will consider ultra-short-term predictions with the length up to several days 
because the method discussed is mostly effective just for such a kind of prediction. Besides, 
ultra-short-term prediction is evidently the most required product. 

The results presented in Fig. [2] clearly show that generally combined predictions provide 
high quality result, more accurate than input series. All three combined series show very 
similar accuracy for the PM predictions. Maybe C3 solution has a marginal advantage. 

However, we have a much more complicated situation with the UT1 predictions. The 
average (CI) prediction does not provide a satisfactory result, probably because of presence 
of bad predictions seen in the central and top- left part of the plots related to UT1. In this 
case, C3 (best yesterday prediction) is clearly the best for UT1 prediction. 

The second test was performed using the actual ERP predictions made at the U.S. 
Naval Observatory, Russian State EOP Service, and Pulkovo Observatory in the period 
from January 2009 through June 2010. Using these three input prediction series, three 
combined prediction series CI, C2 and C3 were computed as described above, and compared 
with the IERS C04 EOP series. Results of comparison, RMS and maximum differences of 
input and combined series with C04 are depicted in the plots in Fig. [3j Again, one can 
see an advantage of combined series with respect to the input ones. In this case, however, 
the CI solution (average prediction) has shown the best result, especially for the maximum 
prediction error. 

3. CONCLUSIONS 

In this paper, we have studied three strategies of computation of combined predictions of 
Polar Motion and Universal Time. Both a posteriori and actual input predictions were used 
for the experiment. The obtained results have shown that a combined prediction can provide 
robust accurate operational EOP forecast. Combined prediction is shown to be the most 
effective for ultra-short-term prediction for a length up to few days. However, more detailed 
investigations are needed to make a justified choice between tested combination strategies. 

Two possible ways of practical use of combined prediction can be considered separately 
or jointly: 

1. Combination of predictions computed at the same analysis center using different strate- 
gies. 

2. Combination of predictions produced at different analysis centers. 

Moreover, combined predictions made using different combination strategies can be used as 
input for the second-stage combination. 

Indeed, to achieve the most accurate result using a combination procedure, input pre- 
dictions should be computed using the best methods. Using all the input prediction series 
without preliminary assessment of their quality can dilute the accuracy of final product. 
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Figure 3: The RMS and maximum errors of the input and combined predictions: actual 
input predictions. 
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